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To gain insights into the mechanisms of cell fate decision in a noisy environment, the effects of intrinsic and 
extrinsic noises on cell fate are explored at the single cell level. Specifically, we theoretically define the 
impulse of Clnl/2 as an indication of cell fates. The strong dependence between the impulse of Clnl/2 and 
cell fates is exhibited. Based on the simulation results, we illustrate that increasing intrinsic fluctuations 
causes the parallel shift of the separation ratio of WhiSP but that increasing extrinsic fluctuations leads to 
the mixture of different cell fates. Our quantitative study also suggests that the strengths of intrinsic and 
extrinsic noises around an approximate linear model can ensure a high accuracy of cell fate selection. 
Furthermore, this study demonstrates that the selection of cell fates is an entropy-decreasing process. In 
addition, we reveal that cell fates are significantly correlated with the range of entropy decreases. 

Cells from across biological kingdoms are continuously engaged in the process of decision-making. A 
growing number of cell types are being described as capable of decision-making under various circum- 
stances. For example, unicellular organisms make vital decisions to enter various phases of the life cycle to 
adapt to environmental changes. In multicellular organisms, precursor cells mature into specialized cell types 
during development. Therefore, the selection of cell fate in response to both internal and external stimuli is 
essential in a cell's life\ As the simplest eukaryote, budding yeast is often used as a model organism to study the 
molecular mechanisms underlying life processes. 

Recently, a quantitative single- cell analysis of commitment dynamics during the mating-mitosis switch in 
budding yeast has been provided^'^. The commitment points are frequently invoked in the explanation of 
differentiation processes. For the mating-mitosis switch process, the purpose of mating is to fuse two haploid 
cells, which must be restricted to the Gl phase prior to the initiation of DNA replication. The point where a cell 
loses mating competence and commits to the cell-cycle is called the Start point (Fig. 1 A). Depending on the 
progression level, the process can be divided into the following two stages: a pre-Start state and a post-Start state. 
As shown in Fig. 1 B, the dynamics can be roughly visualized as a "quasi-potential landscape" in which each 
potential well represents a state. In this pre-Start stage, the process is sometimes reversible to the post- Start state if 
appropriately treated, thus implying instability of the pre-Start statel However, the post-Start state generally 
becomes irreversible if the system passes the critical Start point. As reported previously, the Start point is 
accurately predicted by the nuclear Whi5 concentration independent of cell size, type and Gl duration^. Our 
recent research has also confirmed that Start can be characterized by entropy^, as entropy defines the height in the 
landscape of cell fate decision-making. 

Stochastic fluctuations are ubiquitous in many real dynamical systems as follows: physical, chemical, and 
biological systems. The information transduced in cellular signaling pathways is significantly limited by noise^"^. 
Noise is increasingly appreciated as a force shaping biology. Therefore, the importance of precisely understanding 
the mechanisms of cell fate decision in a noisy environment has already been recognized. Cells may exploit noise 
in different beneficial ways. For instance, phenotypic variability may be triggered by noise because fluctuations 
enable the exploration of the phase space through different types of dynamics^"^\ This variability has been 
observed in several natural systems such as the galactose utilization network in budding yeast^^, the process of 
DNA uptake from the environment in B. subtilis^^'^^, photoreceptor differentiation in the fruit fly retina^^ and 
stem cell differentiation^^"^^. By combining mathematical modeling and gene expression studies in zebrafish, 
Zhang et al. demonstrated that noisy expression can actually facilitate boundary sharpening^^. In addition. 
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Figure 1 | Background of cell fate selection in budding yeast. (A) Schematic of the cell-cycle. (B) Cell fate decision and potential landscape. The pre-Start 
state is situated immediately before the tipping point Start. At this stage, the system is sensitive to external stimuli. The post- Start state is the other stable 
state or a minimum of the potential function where the system is usually irreversible to the pre- Start state. (C) Schematic of the core interacting network 
between the Gl cell-cycle and the pheromone-induced MAPK pathway. All components and reactions considered in our mathematical model are 
included. All arrows for individual reactions are marked with the corresponding number of this reaction in the model. 



interactions among different types of noise have been studied the- 
oretically^"'^^. For example, the advantage of noise in a regulatory 
network in a noisy environment has been demonstrated by Chou et 
^po,2i Interestingly, a general principle states that the capability of a 
feedback system to attenuate the input noise is dependent on the 
difference between the deactivation time and the activation time^^. 

However, due to complexity in highly interconnected biochemical 
networks, many related questions need to be further explored. For 
example, stochastic noise may interfere with the molecular regula- 
tions that cooperatively allow a single cell to choose between two 
different fates. However, the exact roles of stochastic noise on the cell 
fate decision process are not well established. The molecular basis for 
a yeast cell to use such a fundamentally stochastic process to generate 
extremely reliable outcomes also remains unclear. Thus, it is crucial 
to thoroughly characterize the features of stochastic cell fate decision. 
Furthermore, the combined influences of intrinsic and extrinsic 
noises on cell fate decision between cell -cycle commitment and mat- 
ing arrest are poorly understood. 

To select a mitotic fate, the upstream Gl cyclin (cell-cycle signal) 
activates the cyclin- dependent kinase (GDK) to phosphorylate the 
transcriptional inhibitor, Whi5^'^. This initiates a positive feedback 
loop centered on the downstream Gl cyclins, Clnl and Cln2, which 
drives entry into the cell-cycle^^'^^. Conversely, mating arrest is affec- 
ted by Farl, most likely via stoichiometric inhibition of the Gl cyclin 
CDK complexes^'^. Upon pheromone stimulation, Farl is regulated 
in a feedforward manner by the mitogen -activated protein kinase 
(MAPK), Fus3, which interacts physically with the scaffold protein, 
Ste5^^"^^. In addition, the Gl cyclins inhibit the mating pathway by 



promoting the phosphorylation and degradation of both Farl and 
the Ste5 scaffold^'^^. Thus, in this study, we simplified the regulation 
network of the cell-cycle and pheromone-induced MAPK pathway 
as well as the mathematical model presented in the literatures^'^. 
Only the proteins mentioned above, tightly related to the crosstalk 
between the cell- cycle and MAPK pathway, were selected for our 
core interacting network. The schematic diagram used to build the 
mathematical model is shown in Fig. 1 C. We developed a simpli- 
fied deterministic model (Ordinary- Differential-Equation (ODE) 
model) and a corresponding stochastic model (Chemical-Langevin- 
Equation (CLE) modeP^) responsible for selecting cell fates between 
cell-cycle commitment and mating arrest. Based on these two 
classes of models and quantitative measures of cell fate selections, 
we investigated the effects of both external and internal noises on 
the accuracy of cell fate decisions. These findings will provide new 
insights into the mechanisms of cell fate decision in a noisy 
environment. 

Results 

The consistency between the simulation results of the mathema- 
tical models and the experimental observations under wild type 
conditions. To validate our simplified ODE model, the signaling 
dynamics of the system, including the cell- cycle and the 
pheromone-induced MAPK pathways, were simulated. The initial 
values listed in Table 1 were suitable for the cells that had already 
finished at least one cell-cycle, and the cells were all at the beginning 
of Gl phase. Some biomarkers were chosen to indicate different cell 
fates as follows: cell-cycle commitment was represented by high 
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Table 1 | Parameters and initial concentrations of all components 
in the simplified mathematical model 



Parameters 


Value 


Parameters 


Value 


k, 


0.6500 


k. 


10.0000 


k2 


8.0000 


kio 


0.1200 


k3 


0.1000 


k,, 


2.6750 


k4 


0.0310 


ki2 


1 .0000 




0.2600 


k,3 


0.1400 


k6 


0.5800 


ki4 


0.1530 


kz 


0.6200 


di 


0.3000 


ks 


0.3150 


d2 


0.0250 


Signals 


Value 


Signals 


Value 


(X(t) 


2.400(fort>t-added) 


Qi 


0.2500 




or 0 (otherwise) 






Components 


Concentrations 


Components 


Concentrations 




(in the model) 




(in the model) 


Ste5 


1 .0000 


Whi5 


1 .0000 


Farl inact 


1 .0000 







The concentrations of other components that are not listed in the Table are 0. 



levels of Clnl/2 and WhiSP, and mating arrest was predicted by high 
levels of SteSmem and Farl. The detailed comparisons between the 
original model and the simplified model are shown in Supplemen- 
tary text. The numerical results for the simplified ODE model are 
presented in Supplementary Fig. SI. Given the small size of a yeast 
cell (—30 fL) and the low concentration of regulatory proteins 




0 10 20 30 

t[min] 



(—50 nM), the total number of molecules of each regulatory 
protein in a cell is limited (—1,000), and the intrinsic molecular 
fluctuations are not only inevitable but also large enough^^. The 
system volume y was set to be 1000. Thus, the intrinsic molecular 
fluctuation was assigned to a Gaussian distribution with a mean of 0 

and a standard deviation oi \ j Vv ~ 0.036. Fig. 2 A shows the 

simulation result of the CLE model compared to the experimental 
data. Fig. 2 B and C illustrate that the cell would commit without 
pheromone (Fig. 2 B) and arrest with the addition of pheromone 
(Fig. 2 C). The relationship between the activation of WhiSP (i.e., the 
mean proportion of Whi5 exported from the nucleus from 29.8 to 
30 min) and the addition time of pheromone is shown in Fig. 2 D. 
According to our theoretical calculations, the critical ratio of WhiSP 
for WT simulated by our CLE model was approximately 50.43%. 
This result was close to the reported experimental results indicat- 
ing that cell fate is determined when approximately half of the Whi5- 
GFP (52% ± 3%) has been exported^ 

Cell-cycle reentry as simulated by our models. A recent experiment 
confirmed that for cell- cycle reentry from mating arrest, 64% ± 4% 
of nuclear Whi5 is required to be removed compared with 52% ± 3% 
for a cycling cell to progress through Starf. Here, we used a Clnl/2 
signal pulse {sp) (Fig. 3 A) to simulate the exogenous Cln2 pulse used 
in biological experiments. In addition, the degradation rate of Farl^ct 
was set to 0 between 30 min and the end time of the exogenous Cln2 
pulse to maintain a high Farl expression level for the mating cells. 
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Figure 2 | Simulation results of the CLE model. (A) Time course of Whi5P in our CLE model. The experimental data are marked with a star. (B) The time 
course of a WT yeast cell when no pheromone is added. (C) The time course of a WT yeast cell when pheromone is added at 0 min. (D) Relationship 
between the addition time of pheromone and the activation of WhiSP. Here we set R = 0, indicating that the intrinsic fluctuations in all the reaction 
channels are independent of each other. 
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Figure 3 | Simulations of reentry and reversibility in cell fate selection. 

(A) Time course of Clnl/2 signal pulse (sp). (B) The time course of a WT 
yeast cell when pheromone is added at 0 min in the ODE model. The 
duration is 3 min, and the cell fate remains in the mating program. (C) The 
time course of a WT yeast cell when pheromone is added at 0 min in the 
ODE model. The duration is 30 min, and the cell reenters the cell-cycle 
from the mating program. (D) The time course of a WT yeast cell when 
pheromone is added at 0 min in the CLE model. The duration is 3 min, 
and the cell remains in the mating program. (E) The time course of a WT 
yeast cell when pheromone is added at 0 min in the CLE model. The 
duration is 30 min, and the cell reenters the cell-cycle from the mating 
program. 

Similar to Doncic's experimental design^, we varied the durations 
between 3 and 30 min to express variable amounts of exogenous 
Cln2 in our mathematical model. The simulation results are 
presented in Fig. 3 B-D. The pre-Start state cells can clearly 
commit reversibly to the post- S tort state if appropriately treated. 
The relationship between the duration time of the Clnl/2 signal 
pulse and the maximum WhiSP value approximately 10 min after 
the end of the pulse (corresponding approximately to the time 
(5 min) it takes to inactivate the MET3 promoter^°) is shown in 
Supplementary Fig. S2. The critical point that had the largest 
derivative in the dose -response curve had a critical WhiSP amount 
of 65.42%. This result agreed with the experimental result of 64% ± 
4%. 

Cell fates can be distinguished by the impulse of Clnl/2. Based on 
the experimental data, the histogram of WT yeast cells is shown in 
Supplementary Fig. S3. The histogram (the number of cells vs. the 
fraction of exported Whi5) indicates the number distribution of cells 
selecting different fates. The bin size of our histogram was set to 0.1. 
As mentioned in Ref. (2), the nuclear Whi5-GFP concentration was 
sufficient to predict cell fate in approximately 97% of WT Gl cells. 
Supplementary Fig. S3 shows that the number of cells possessing an 
arrest fate slightly increased as the fraction decreased from 
approximately 0.6 to 0.1 and that the number of cells committing 
increased as the fraction increased from approximately 0.6 to 1. It 
should be noted that this property was not captured by most of our 
numerical simulations. In the experiment, the cell fates of 
approximately 280 cells were observed, and the theoretical analysis 
was performed using the data of these 280 cells. Due to the lack of 
more experimental details, the concentrations of each protein in 
different cells were unclear. In addition, the study was focused on 
the coordination of intrinsic and extrinsic noises in cell fate 



decisions, and the initial concentration distribution would not 
change the roles that the noises play. Therefore, in our theoretical 
model, we considered the initial concentration of WhiSP from the 
different cells as a uniform distribution without loss of generality. 
However, this interesting property can also be reproduced by our 
simulation result if the initial concentrations of WhiSP for different 
cells are sampled from the probability distribution of concentrations 
in Fig. S3. The simulation result is shown in Supplementary Fig. S4. 

The histograms of the simulation results of the ODE model (Fig. 4 
A and B) and CLE model (Fig. 4 C and D) are shown. Each cell is 
denoted by different random initial values in the pheromone envir- 
onment. We utilized two numerical indexes to distinguish different 
cell fates in budding yeast, (i) In Fig. 4 A and C, we used the fraction of 
WhiSP after the Gl duration, i.e., the average fraction of WhiSP at 
29.8-30 min, to identify different cell fates. If the amount of WhiSP 
from our model was larger than the critical amount of WhiSP at 
Start, the cell would commit; otherwise, the cell would arrest, (ii) 
In Fig. 4 B and D, the other numerical index was the impulse of Clnl/ 
2 (See Methods) where the cell would commit if the impulse was 
larger than the critical impulse otherwise, the cell would arrest. The 
comparison between Fig. 4 A and Fig. 4 C (or Fig. 4 B and Fig. 4 D) 
revealed that the WhiSP index was inadequate to the intrinsic noise 
but that the impulse index benefitted from the fluctuations. Hence, 
under the perturbation of noise, the impulse of Clnl/2 can provide a 
reliable discrimination for cell fates due to the time integral of Clnl/2. 
In addition, the impulse of Clnl/2 enabled a more effective separa- 
tion, with an accuracy of 99.17% (Fig. 4 D) compared with the frac- 
tion of WhiSP after the Gl duration. The accuracy of 99.17% was 
slightly better than the value of approximately 97% derived from the 
experimental data in Ref. (2). The histogram in Fig. 4 D conformed to 
the experiment results shown in Supplementary Fig. S3. The clas- 
sification result validated our assumption that cell fates can be dis- 
tinguished by the impulse of Clnl/2. This phenomenon suggests that 
the continuous operations of Gl cyclins, but not the instantaneous 
concentrations, determine the cell fates of yeast and clearly indicates 
that the impulse helps to decrease the uncertainty of cell fate at a 
molecular level. Therefore, the critical impulse of Clnl/2 was selected 
as the discrimination factor between mating and cell-cycle in the 
following sections. The results of R = 0.2S and R = 0.7S are shown 
in Supplementary Figs. SS and S6. The result of R = 0.2S was similar 
to the uncorrected result, but the higher correlation coefficient 
caused a competitive cell-cycle process. These results implied that 
the yeast cells adapted to the lower correlation between internal 
noises, which corroborated the logic of intrinsic noise generally being 
specific to any one gene. Therefore, the correlation coefficient of 
intrinsic noise R was set to 0 in the following sections. 

Large extrinsic noises and uncorrelated extrinsic noises disturb 
the separation of cell fates. Extrinsic noise is common to genes of 
any one cell due to differences between cells, such as temperature and 
energy state^\ When extrinsic noise is added, each parameter in the 
model is sampled from a Gaussian distribution with a mean equal to 
the listed value and a standard deviation equal to a. The following 
three types of extrinsic noise were considered: (1) the extrinsic signal 
noise of the mating pathway, i.e., the extrinsic noise of pheromone 
introduced in the parameter of a(t) where the noise intensity of the 
signal in the mating pathway is denoted as abating; (2) the extrinsic 
signal noise of the cell- cycle pathway considered in the parameter of 
ai, where the corresponding noise intensity is denoted with acydeJ 
and (3) the extrinsic noise of other biochemical reaction rates, where 
the noise intensity is marked with Qbrc)- And we would use 'a' to label 
the standard deviation of all types of extrinsic noise without causing 
ambiguity. We also considered the correlation between each extrinsic 
fluctuation. Here two cases were adopted to simplify the situation: 
correlation coefficient Rgx is equal to 1 or 0. In Fig. S, the correlation 
coefficient Rex was set to 1, and the noise intensity a was set to 1 (the 



SCIENTIFIC REPORT: | 4 : 5764 | DOI: 1 0.1 038/srep05764 



4 



40 
35 
30 
25 
z 20 
15 
10 
5 
0 



I arrest 
H commit 



0.1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 
Fraction of Whi5 removed from nuclear 



40 

35 
30 
25 
20 
15 
10 
5 



I arrest 
H commit 



0.1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 
Fraction of Whi5 removed from nuclear 



40 
35 
30 
25 
20 
15 
10 



I arrest 
H commit 



40 
35 
30 
25 
20 
15 
10 



I arrest 
H commit 



JL 



0.1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 
Fraction of Whi5 removed from nuclear 



0.1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 
Fraction of Whi5 removed from nuclear 



Figure 4 | Simulation histograms based on cell fate determined by the fraction of exported Whi5 at the time of pheromone addition. N indicates the 
number of cells. (A) The cell fate is distinguished by Whi5P after Gl duration in the ODE model. (B) The cell fate is distinguished by the impulse of Clnl/2 
in the ODE model. (C) The cell fate is distinguished by Whi5P after the Gl duration in the CLE model. (D) The cell fate is distinguished by the impulse of 
Clnl/2 in the CLE model showing a perfect threshold system. We set R = 0 for (C) and (D). 



results of extrinsic the noise intensity of 0.1, 0.2 and 0.5 with Rex = 1 
are shown in Supplementary Figs. S7-S9, respectively). When 
comparing Fig. 5 A, to Fig. 4 B, we found that the large extrinsic 
noise of the parameters would promote the competitiveness of 
mating pathway but fail in the separation of cell fates. If the large 
extrinsic signal noise of mating pathway was also taken into 
consideration, we found that the extrinsic noise of pheromone 
would promote the competitiveness of the mating program and 
disturb the separation of cell fates (Fig. 5 B). Fig. 5 C indicates that 
the large extrinsic noise of the cell-cycle pathway would enhance the 
competitiveness of the cell- cycle pathway but would not help the 
separation of cell fates. The result of all the extrinsic noises are 
shown in Fig. 5 D, it suggested that the large extrinsic signal noise 
of the cell- cycle pathway is more aggressive compared to the signal 
noise of the mating pathway with the same standard deviation. 

The results of the correlation coefficient R^x ^ 0 and extrinsic 
noise intensity a = 1 are illustrated in Supplementary Fig. SIO 
(the results of extrinsic noise intensity of 0.1, 0.2 and 0.5 with Rgx 
= 0 are shown in Supplementary Figs. SI 1 -SI 3). Supplementary Fig. 
SIO to Supplementary Fig. SI 3 suggested that uncorrected extrinsic 
noises cause more misclassification and more cross-borders in dif- 
ferent cell fates. These results suggested that cells can adapt to highly 
correlate extrinsic noises, and these results corroborated our general 



viewpoint that the origin of extrinsic noises, such as temperature, has 
a similar influence on each reaction channel. Therefore, the uncor- 
related extrinsic noises will disrupt the cell fate decision due to the 
lack of adaptation to these noises. The correlation coefficient Rex was 
set to 1 in the following sections. 

We then considered that the extrinsic signal noises in the cell -cycle 
and mating pathways have different strengths compared with noise 
in other biochemical reaction parameters. The results of the signal 
noise with intensities of 0.1 and 1 in addition to the noise in each 
reaction parameter with an intensity of 0.5 are shown in 
Supplementary Figs. S14 and SI 5. Figs. S9, S14 and SI 5 show that 
more cells became committed with greater strengths of extrinsic 
signal noise in the mating pathway and that the extrinsic signal noises 
played a more important role than the extrinsic noise of the reaction 
parameters. This property indicated that the extrinsic signal noises 
would greatly confuse the fate selection, and this confusion may 
result in the trigger role that the signals play. 

In addition, we explored the role of the extrinsic signal noise without 
the extrinsic and intrinsic noises of the reaction parameters (detailed in 
Supplementary Figs. S16 and SI 7). Compared with each of the noises 
alone, we concluded that a small fluctuation in pheromone together 
with a relatively large fluctuation in the cell-cycle signal give rise to a 
higher performance in promoting the accuracy of fate selection. 
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Figure 5 | Simulation histograms based on cell fate as determined by the fraction of exported Whi5 at the time of pheromone addition. The cell fate is 
distinguished by the impulse of Clnl/2. (A) The result of the ODE model with extrinsic noise in the biochemical reaction parameters (abrc = !)• 

(B) The result of the ODE model with extrinsic noise in the biochemical reaction parameters and the mating pathway signal (abrc = 1 and amating = !)• 

(C) The result of the ODE model with extrinsic noise in the biochemical reaction parameters and signal of the cell-cycle pathway (abrc — 1 and acycie = 1 )• 

(D) The result of the ODE model with extrinsic noise in the biochemical reaction parameters, the signal of the mating pathway and the signal of the cell- 
cycle pathway (abrc = 1> CTmating = 1 and acycie =1). Rex = 1 indicates that the correlation coefficient of the considered types of extrinsic noise equals 1. 



Intrinsic noise and extrinsic noise play different roles in cell-fate 
selection. In this experiment, the intrinsic and extrinsic noises in 
yeast cell populations were integrated. At first, the system volume 
was set to 1000, and the effects of extrinsic signal noise were then 
examined. The results of the cell-cycle signal noise and the mating 
signal noise are shown in Supplementary Figs. SI 8 and SI 9. In 
contrast to the results shown in Fig. 4 D, the cell-fate selection 
depends almost entirely on intrinsic noise for a small extrinsic 
noise with a standard deviation of 0.1 and 0.2 (Supplementary Fig. 

518 A and B; Supplementary Fig. S19 A and B). Supplementary Fig. 
SIS C and D as well as Supplementary Fig. S19 C and D illustrate that 
the selection of cell-fate was influenced by strong extrinsic noise. 
Furthermore, strong extrinsic noise in the mating signal caused a 
large decrease in the accuracy of fate decision (Supplementary Fig. 

519 D). Previous studies have suggested that noise is dominated by 
intrinsic fluctuations at low expression levels. Thus, based on this 
conclusion, the results shown in Supplementary Figs. SI 8 and S19 
confirmed that the extrinsic noise values were close to a well- 
delineated CVext less than 0.5, which was in agreement with 
previous studies^^"^^. In addition, the effects of intrinsic noise were 
studied for different system sizes (V = 10000, V = 1000, V = 100, 



and V = 10) with a fixed standard deviation of extrinsic noise. The 
corresponding histogram results with the fixed standard deviations 
of 0.1 and 1 for the extrinsic noise are plotted in Fig. 6 and 
Supplementary Fig. S20, respectively. Based on the observations 
above, increasing intrinsic fluctuations resulted in a parallel shift in 
the separation ratio of WhiSP, and the increasing extrinsic 
fluctuations led to the mixture of different cell fates. This result 
was consistent with the logic of extrinsic fluctuations being an 
unwelcome influencing factor while cells adopt the intrinsic 
fluctuations to decease the confusion in cell fates. 

Extrinsic noise in Clnl/2-mediated phosphorylation of Whi5 
plays an important role in the separation of cell fates. In this 
section, we focused our attention on the function of a single 
extrinsic or intrinsic noise. As mentioned above, extrinsic and 
intrinsic noises play different roles in the selection of cell fate, but 
it remains unknown which one (or ones) acts as the most important 
factor. As shown in Fig. 7, the extrinsic noise in ki2 with both small 
and slightly large standard deviations benefitted the mating process 
(Fig. 7 A and B), but the noise with a larger standard deviation led to 
an orderless separation of cell fate (Fig. 7 C and D). However, 
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Figure 6 | Simulation histograms based on cell fate as determined by the fraction of exported Whi5 at the time of pheromone addition. The cell fate is 
distinguished by the impulse of Clnl/2, and extrinsic noise has a standard deviation a = 0.1. The intrinsic noise has a standard deviation (A) 1 / W ~ 
0.01 (V = 10000). (B) 1 /W 0.0316 (V = 1000). (C) I / Vv ^ 0.1 (V = 100). (D) 1 /v^ « 0.3162 (V = 10). 



extrinsic noises of other parameters (taking kio as an example) will 
play a similar role even though the standard deviation was quite large 
as shown in Supplementary Fig. S21. Hence, our simulation results 
illustrated that only the extrinsic noise from ki2 (a Michaelis 
constant) induced the mixture of different cell fates indicating that 
the extrinsic noise from the reaction of Clnl/2 -mediated 
phosphorylation of Whi5 plays an important role in the separation 
of cell fates. This result suggested that a potential mechanism in 
decreasing the extrinsic noise of ki2 is essential for budding yeast. 
In fact, this corresponding reaction of ki2 (i.e., Clnl/2 -mediated 
phosphorylation of Whi5) only occurs in the nucleus to reduce the 
extrinsic noise. The roles of single intrinsic noises in each reaction 
channel are presented in Supplementary Fig.S22. By numerical 
analysis, we classified all the reaction channels into two categories 
as follows: supporter of the mating pathway (resulting in more 
mating cells than Fig. 4 B) or cell -cycle program (resulting in more 
committed cells than Fig. 4 B). Interestingly, interaction between 
Clnl/2 and Farl^ct acted as the most important balance in the cell 
fate decision. It is clear that the intrinsic noise in the inhibition of 
Clnl/2 on Far 1 act belonged to the first three supporters of the mating 
program, and the inhibition of Far 1 act on Clnl/2 was the unique 
supporter of the cell-cycle program. It is well known that the 
interaction between Clnl/2 and Far 1 act is a crucial crosstalk 
between the cell -cycle and MAPK pathway. Our findings further 



verified that the noise from the crosstalk also contributes to the 
fate decision. 

The linear interplay of intrinsic and extrinsic noises ensures a high 
accuracy of cell fate selection in budding yeast. Furthermore, to 
quantitatively clarify the relationship between extrinsic and intrinsic 
noises, we proposed a formula to quantify the accuracy of cell fate 
selection (See Methods) and calculated Acc under different noises. 
The numerical results are shown in Fig. 8. Fig. 8 A shows that for 
small intrinsic noise, the accuracy of cell fate selections increased at 
first and then decreased with the increase of standard deviation of 
extrinsic noise. The result remained unchanged for small extrinsic 
noise when increasing the standard deviation of intrinsic noise. We 
also identified the point(s) with a maximum accuracy larger than 
95% for each fixed intrinsic noise. The first three largest extrinsic 
noise strengths among these points for each fixed intrinsic noise were 
taken into consideration (only the largest one will be adopted if the 
total number of points is less than three). The average values of these 
three extrinsic noise strengths are marked with stars for each row in 
Fig. 8 A. Based on our analysis, these points were found to be 
approximated by the following straight line: yi = —1.2391 — 
0.1235X. The linear regression showed that it had the property of 
0.8076 in the Pearson correlation coefficient suggesting that a proper 
strength of extrinsic and intrinsic noise near the straight lineyi helps 
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Figure 7 | Simulation histograms based on cell fate as determined by the fraction of exported Whi5 at the time of pheromone addition. The cell fate is 
distinguished by the impulse of Clnl/2, and only extrinsic noise in ki2 is taken into consideration. (A) a = 0.05. (B) a = 0.1. (C) a = 0.5. (D) a = 1. 
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Figure 8 | Selection accuracy in noisy environment. (A) The accuracy as functions of the strength of extrinsic noise a and intrinsic noise 1 /W. The 

average values of the three largest extrinsic noise strengths are marked with stars. The straight line yi is obtained using linear regression. (B) The contour 
lines of the accuracy function shown in (A). The points with stars represent the near edge points in each contour Hne, and these data points are fitted well 
by a straight line y2. 
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Figure 9 | Cell fate and entropy. (A) The fraction of entropy- decreasing cells for mating and cell-cycle fates. (B) The fraction of cells with an entropy 
fluctuation over or below 15% for mating and cell-cycle fates. 



to ensure a high accuracy in budding yeast cells. The contour lines of 
Fig. 8 A are shown in Fig. 8 B. The contour lines had the approximate 
right angle form. The stars in Fig. 8 B indicated the near edge points 
in each contour line, and these points can also be approximated by 
the following straight line: y2 = 0.7288 + 1.0769x. Therefore, the 
accuracy more likely depended on the intrinsic noise for the tuple 
(logjQ(l/v^),logjQ((7)) under the line and the accuracy was 
mainly determined by the extrinsic noise for the upper tuple 
(logjo(l/v^),logjo(o-)). In addition, the cross point of these two 
straight lines was ( — 1.6394,-1.0366), which was close to the 
default noise strength of (logio(l /\/lOOO),logio(0. 1)) = ( - 1 .5, - 1) 
obtained through the experiments. These results indicated that the 
noise strengths in wild type yeast cells contribute to ensure a high 
accuracy in fate selection as an inevitable consequence of intelligent 
evolution. 

Fluctuation range in entropy shows strong correlation with the 
cell fates. Network entropy exhibits dynamic changes in time course 
differentiation data, and it is in line with the cell fate stage^^. We 
recorded the entropy of a yeast cell Y between 0 and 30 min every 
0.01 min, and the entropy was represented by Ey. The data were 
collected through the CLE model, and we set a = 0.1 and V = 
1000. For a yeast cell Y, we further computed the ratio of an 
entropy- decreasing time point, i.e., the time point where entropy is 
larger than the entropy of its following neighbor. The result is shown 
in Fig. 9, and it suggested that approximately 92.14% of the cells had 
more than 1500 entropy- decreasing time points. These results 
suggested that the cell fate selection program may be an entropy- 
decreasing process. Thus, we investigated if there were any 
differences in the fluctuation range between different cell fates. The 
fluctuation range for a yeast cell Y was defined as follows: 

max(£Y)-niin(EY) , . r n 

Fy = , where Eq is the entropy of cell Y at 

0 min. As observed in Fig. 9 B, the fluctuation range showed 
strong correlation with the cell fates (R = 0.8522 and p-value < 
10~™). This result illustrated that the entropy of these cells 



selecting the mating fate was decreased significantly due to the 
input of pheromone information. However, these cells selecting 
cell -cycle fate shutdown the input of the mating signal and 
underwent a small fluctuation in entropy. This result was 
consistent with the theory of information entropy and negative 
entropy principle. In another way, the result suggests that the cell 
fate is determined by the impulse of Clnl/2 is somehow reasonable. 

The correlations between entropy-decrease and cell fate for the 
CLE model with other tuples (logjQ(l/V^),logjQ((7)) characterized 
by the two linear relationships in Fig. 8 are shown in Supplementary 
Fig. S23. These results also suggested a great dependence between cell 
fate and the decrease of entropy, thereby indicating that a decrease in 
entropy can serve as a common rule for the selection of cell fate in 
budding yeast. In addition, the dependence of the average fluctuation 
in entropy for right predicted cells and the noise characterizing 
model was calculated and shown in Fig. 10. Fig. 10 shows that the 
accuracy decreased with an increase of logjQ(l / v^) according to the 
linear model y2. The high accuracy tuples (logjQ(l/A/y),logjQ((7)) 
characterized by model yi, model y2 as well as the combinatorial 
model of yi and y2 (i.e., in the left plane) shared a sharp boundary 
in entropy fluctuation for different cell fates (the left plane) while low 
accuracy tuples suffered an unsharp boundary (the right plane). 
Thus, we conclude that cells in the high accuracy plane have smaller 
fluctuations and consistent differences in entropy decrease while the 
cells with low accuracy in the right plane have larger fluctuations and 
confused differences in entropy decrease, thereby suggesting that 
high accuracy results from both the sensitive decrease in entropy 
and the sharp boundary in entropy fluctuation. 

Sensitivity analysis of the model parameters. In this section, the 
CLE model with a = 0.1 and V = 1000 was used to perform the 
calculations. We perturbed each parameter with 10-30% variations. 
The results of the sensitivity analysis are presented in Table 2. The 
sensitivity analysis implied that ai (the signal of the cell-cycle), ki4 
(Whi5-IClnl/2) and kn (Whi5 + Clnl/2 Whi5P) had relatively 
significant effects on Whi5P expression. The other parameters had a 



SCIENTIFIC REPORT: | 4 : 5764 | DOI: 1 0.1 038/srep05764 



9 



■D 

a> 



-0.05 



-0.1 



S -0.15 



Q. 
O 



(U 

> 

< 



-0.2 



5= -0.25 



-0.3 



-3 



The left plane 



C-correct && Model 

C-correct && Model 

C-correct && Model (y^+y2)/2 

—I M-correct && Model y^ 

-0 M-correct && Model y^ 



— M-correct && Model (y^+y2)/2 




-2.5 



-1.5 



-0.5 



Figure 10 | The dependence of entropy decreasing and accuracy. Here, C-correct indicates cells with more than 50% of the Whi5 removed from the 
nucleus, which are predicted to commit, and M-correct indicates cells with more than 50% of the Whi5 still in the nucleus, which are predicted to arrest. 



slight effect or almost no effect on WhiSP expression. These results 
demonstrated that the Start point was mainly dominated by the 
biochemical processes of Whi5. In addition, our model showed 
good dynamical robustness. Furthermore, we perturbed all of the 
parameters for 10%, 20% and 30%, and the simulated results 
indicated that the number of cells predicted to commit varied 
slightly with values of 12.24%, 12.59% and 7.14%, respectively. 

Discussion 

Despite considerable advances in the research of underlying mechan- 
isms in cell fate decision, it has been largely unclear how individual 
cells make critical cell fate decisions in the presence of stochastic 
fluctuations in intracellular biochemical reactions and extracellular 
environments. The relative importance of intrinsic and extrinsic 
noises in cell fate decision also remains uncertain. We attempted 
to answer these questions with a simple theoretical model of the cell 
fate decision network responsible for cell-cycle commitment and 
mating arrest in budding yeast. Using ODE and CLE models, the 
combined effects of extrinsic and intrinsic noise on cell fate selection 
were clarified at a single cell level. As a driving factor of the Start 
point, the impulse of Clnl/2, which can effectively discriminate 
between different cell fates, was used to identify cell fates. In this 
study, we found that the uncorr elated extrinsic noises caused mis- 
classification and cross-borders in different cell fates. Increasing 
intrinsic fluctuations gave rise to a parallel shift in the separation 
ratio of Whi5P, and increasing extrinsic fluctuations resulted in a 
mixture of different cell fates. Our investigation of the function of 



single extrinsic noise or intrinsic noise suggested that an effective 
method to decrease the extrinsic noise in the Michaelis constant of 
Clnl/2 -mediated phosphorylation of Whi5 is required for budding 
yeast. Furthermore, our study indicated that the linear interplay of 
intrinsic and extrinsic noises contributes to a high accuracy of cell 
fate selection. It should be pointed out that the biochemical reaction 
network is highly nonlinear. Our key findings, i.e., that the approx- 
imate linear coordination between extrinsic and intrinsic noise may 
help cells transit more accurately from one fate to the other fate, were 
still a result of the nonlinear nature of biological processes. In addi- 
tion, our study also suggested that the selection of cell fates is an 
entropy- decreasing process. We found that a large range of entropy- 
decreasing correlated with mating arrest and that a small range of 
entropy- decreasing was related to the fate of cell- cycle commitment. 

The effects of noise on cell fate decision were quantitatively cla- 
rified by our theoretical model. Conventional scientific thinking has 
been that the random nature of such fluctuations within cells inter- 
feres with the reliable operation of biological systems. Our findings 
showed that the cells could still select the different cell fates due to the 
stochastic noise even if the cells had the same initial conditions. The 
coordination between extrinsic and intrinsic noise was beneficial by 
helping cells transform more accurately from one fate to another. 
Therefore, our work provides a new paradigm suggesting that cel- 
lular noise might have a constructive function in cell fate decision. In 
addition, by controlling the noise level within the yeast cells, we could 
prevent or prompt yeast cells to transition into a specific fate type by 
essentially tuning cellular behavior. 



Table 2 Sensitivity analysis of parameters in the CLE model 


Parameter (Pi) 


Reaction 


APi = 10% 


APi = 20% 


APi = 30% 


Average 


ai 


The signal of Cell- Cycle 


0.8231 


0.8396 


0.8790 


0.8472 


ki4 


Whi5 1 Clnl/2 


0.4450 


0.4476 


0.4512 


0.4479 


kii 


Whi5 + Clnl/2^ Whi5P 


0.2429 


0.2440 


0.2456 


0.2442 


kl3 


Whi5P^ Whi5 


0.0827 


0.0827 


0.0827 


0.0827 


kio 


The degradation rate of Cin 1 /2 


0.0746 


0.0747 


0.0747 


0.0747 


ks 


Self-promoting rate of Cin 1 /2 


0.0100 


0.0100 


0.0100 


0.0100 


others 




<0.01 
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The budding yeast is a relatively well-studied example of a network 
with noisy environments. The cell-to-cell variability of the cell-fate 
decision system and the mating pheromone response pathway in 
yeast have been studied experimentally^^. The quantitative output 
of the cell-fate decision system in response to the defined perturba- 
tions has been estimated. The mechanism that allows cells to respond 
more precisely to pheromones in the presence of a large variation has 
been provided. Furthermore, a combination of genetic elements, 
including pathway regulators, RNA-based transducers, constitutive 
promoters and pathway-responsive promoters, has been used to 
build modular network diverters^^. Cell fate has been engineered 
through synthetic gene circuits, and cell fate determination has been 
measured through halo assays. Motivated by these experiments, it is 
expected that our theoretical results could be tested by natural and 
synthetic genetic circuits in individual cells in the future. 

Our research presents a theoretical framework for cell fate 
decision in budding yeast and provides a clear investigation of the 
interaction between intrinsic and extrinsic noises. Understanding the 
mechanisms that increase the accuracy in cell-fate decision system 
might enable new therapeutic interventions. However, some pro- 
blems should be discussed. First, the dynamics of a biological system 
can be modeled more accurately using a Markov jump process^^"^°. 
Exploring the Markov jump process of this decision-making system 
using the Gillespie algorithm^^"^^ may provide a molecular realiza- 
tion that is much closer to reality. Second, cellular transformations 
that involve a significant phenotypical change of the cell's state often 
use bistable biochemical switches as underlying decision sys- 
tems^°'^^'^^. Hence, the bistability and irreversibility of our simple 
model should be clarified. Third, the relationship between the accu- 
racy of the cell fate and entropy- decreasing process is not studied. 
Our next task will be to analyze the change of entropy in the high 
accuracy region. Fourth, the regulatory mechanism to yield high 
accuracy of cell fate decision in large populations of cells is still poorly 
understood despite the underlying stochasticity^^'^^'^^. The potential 
roles of positive feedback, feedforward, slow transcription and rapid 
phosphorylation in the reUable fate decision needs to be further 
analyzed^'^^. 

Methods 

Mathematical models. To clarify the mechanism for selecting cell fates between cell- 
cycle commitment and mating arrest in a noisy environment, we presented two 
classes of mathematical models. First, based on mass action law and Michaelis- 
Menten kinetics, the dynamics of the regulation network of cell- cycle and 
pheromone-induced MAPK pathway presented in Fig. 1 C were characterized by a 
deterministic model, i.e., the ODE model. For simplification, the activities for total 
Ste5 and total Whi5 were assumed to be one unit. The parameters are listed in Table 1 . 
We estimated these parameter values to reproduce the biological observations in 
budding yeast and the critical transition value for Start in the previous experiment 
simultaneously. To consider the effects of intrinsic noise, the corresponding CLE 
model based on the above ODE model was as follows: 



d[Ste5] 
dt 



- a{t)[SteS] + h [Ste5]^^^ [Clnl/2] + h [Ste5]^^^ - d, [SteS] 

^ V^^OIs^CiCO+^Tms^^ (1) 



)= ^jk^St^^Ut) - ^d,[Ste5]Ut) 



= a{t)[Ste5] - k, [SteS]^^^ [Clnl/2] - d, [SteS]^^^ + ^ V^^Ofs^CiCO 
, ^ (2) 

1 — \ — 7~. — TT. . . 1 ■ 



^h[Ste5]^^janl/2\Ut) - ^ ^ d,[Ste5]^^^Ut) 



dt 



^ ^JHmrnern[F^rl\^actUt) + ^ ^J h[Ste5]^emUt) (3) 



d[Farll 
dt 



- h [Ste5]^em [Farll^act - h [Far 1],,, [Clnl/2] - h [Farl],,, - ^2 [Farl],,, 

- ^ \jHSte5]^em[Fcirl],,acUt) ' ^ ^^6 [Far l],,JClnl/2] ^9(0 (4) 

- ^ ^Jh[Farl]^^Mt) - ^ V^2[Farl],„,,,Cio(0 



d\C\nl/l] , , n r , / n A;8[Clnl/2] 
^ ' ^ --/c7[Farl],,JClnl/2]+ai ' ' ^ 



dt 



^9 + [Clnl/2]^ 

/cio[Clnl/2]-/ci4[W/z/5]--^^/c7[Farl],,JClnl/2]Cn(0 



- V^io[Clnl/2]Ci4(0 - ^ V/ci4[W/2/5]Ci5(0 



(5) 



d[Whi5\ _ /cn[W/z/5] [Clnl/2]^ ^ku[Whi5P] - / ^n[^^^"5] [Clnl/2]^ , ^^^^^ 



dt 



A;i2 + [Clnl/2]' 
+ -^v'^i3[W/z/5P]Ci7(0 



Vvy ^12 + [Clnl/2]' 



d[Whi5P\ _ [W/zz5] [Clnl/2]' 



dt 



^12 + [Clnl/2]' 



-ku[Whi5P]- 



1 



A;n[W/zz5][Clnl/2]^ 
A;i2 + [Clnl/2]' 



(6) 



(0 

(7) 



zV%3[W^Cl7(0 



where is Gaussian white noise, for / = 1, 2, 3,. . .,17, <^i{t)> = 0, and for each /, j 
= 1,2, 3,. ..,17, 



(R§(t-t'),fori^j 
'{ d{t-t')Jori=j ■ 



(8) 



here is a constant number and V denotes the total number of molecules of each 
regulatory protein. Furthermore, when extrinsic noise was considered, each 
parameter pi in the model was sampled from a Gaussian distribution with a mean 
equal to the listed value and a standard deviation equal to a, i.e., pi{l + (JEi{t)), here 
Bi{t) is Gaussian white noise, <£,(t)> = 0, and for each 



ReJ(t-t'),fori^j 
3{t-t'),for i=j 



(9) 



where Rex = 1 indicates that the extrinsic noises we considered were correlated with 
each other; and Rgx = 0 indicates that all these noises were independent. 

The CLE model was solved with the Euler-Maruyama method"^"^, and the temporal 
step size was set to 0.01 min. We simulated different yeast cells by stochastic initial 
values. Each cell had a set of random initial values. We obtained the initial concen- 
tration of WhiSP randomly at first, and the initial concentrations of the remained 
components are given according to the interactions in the simplified model. A. 
Doncic et al. compared the amount of Whi5-GFP exported to the largest value of 
nuclear Whi5-GFP within 30 min prior to pheromone exposure for each Gl cell. We 
recorded the amount of Whi5 exported, i.e., the stochastic initial value of WhiSP. 
Thus, we simulated the histograms based on cell fate through a theoretically defined 
cell fate. 

Determination of the Start point. In yeast, the Gl checkpoint Start determines 
whether a cell enters the mitotic cycle or engages the mating program but never both 
because the aims of the two programs are diametrically opposed as follows: mating 
produces one cell from two; and mitosis produces two cells from one. Doncic et al. 
identified the main probe for the status of the Start trigger as the amount of WhiS 
protein exported from the nucleus^. Our previous study showed the relationship 
between the activation of WhiSP (i.e., the proportion of WhiS exported from the 
nucleus at 30 min calculated by our full ODE model) and the addition time of 
pheromone, and a critical point that had the largest derivative in the dose-response 
curve was defined to be the Start point. The Start point defined above was proven to be 
reasonable. 

Impulse of cell fate. Because the cell fate can only be determined experimentally, the 
driving factor of Start must be defined if we want to trace back to the underlying 
mechanism. It is known that the Gl cyclins collaborate to drive yeast cells through the 
Gl-S transition. According to classical mechanics, the impulse is the multiplication of 
the applied force and time interval equal to the change in momentum. Hence, the 
impulse will cause a change in mechanical motion of the particle. The impulse of 
Clnl/2, which was defined as the accumulation of Clnl/2 with time, was presented as 
follows: 
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lit .added) = 



Clnl /2t_added(t)dt 



The relationship between the Clnl/2 impulse and the addition time of pheromone 
was provided. The impulse of the point that had the largest derivative in addition time 
v.s. the proportion of exported Whi5 curve is further defined to be the critical impulse 
of cell fate. When a cell treated with the a-factor at 0 min was taken into consideration 
in our ODE or CLE model, the impulse of this cell was determined. We verified the 
assumption that the cell would commit if its impulse was larger than the critical 
impulse or otherwise arrest. 

Entropy defined in our network. To further quantify the dynamics of the network, 
we provided a definition of the entropy E in this biochemical reaction network by 
extending the concept of entropy in statistical mechanics as follows: 



(11) 



£=-^P(01og(P(0) 

i=l 

where P(i) = x(t)/( ^ x(j)) is the distribution function of the proteins' 

;=1 

concentration; x{i) denotes the concentration of the zth component, and n is the 
number of all components considered in the model (Here n = 7). From the viewpoint 
of evolutionary biology, the change of the entropy implies the evolution of the 
structure and function of a molecular network due to its adaptation ability. 

Quantification of cell fate's accuracy. We observed that if the mating pheromone 
was applied to a cell when more than 50% of Whi5 was still in the nucleus, the cell 
almost always directly arrested, which is called pre-Start. If the mating pheromone 
was applied when more than 50% of nuclear Whi5 was removed, the cell almost 
always underwent one more mitotic cycle before arresting, which is called post-Start. 
Thus, the accuracy Acc can be defined as follows: 



Acc = 



Nm 



Total 



(12) 



where NM-correct denotes the number of cells with more than 50% of the Whi5 still in 
the nucleus, which are predicted to arrest by our critical index; and Nc-correct denotes 
the number of cells with more than 50% of the Whi5 removed from the nucleus, 
which are predicted to commit. The total number of cells in our model was set to 280, 
which was consistent with the average number of cells used by Doncic et aP. 

Measure of sensitivity and robustness. Due to the lack of experimental data to 
determine all of the parameters, it was necessary to analyze the sensitivity of the 
system to changes of the parameters. Similar to the literature (4), each parameter was 
perturbed with 10-30% variations. The sensitivity function Sj{t) of parameter Pj at 
time t was defined as follows: 



1 0(Pj + AP.- , 0 - 0{Pj - APj , t) I 

dOit) dPj{t) J oiP„t) 

0(0 / Pjit) 



where 0(t) is the model output (Whi5P level) at time t; Tis the total reaction time; AP^ 

is a small perturbation; and Sj = Sj{t)dt is the sensitivity value of parameter Pj. In 
Jo 

addition, we also perturbed all the parameters with 10-30% variations and measured 



the robustness as follows: 



RM(AP) = 



\C{P + ^P)-C{P-^P)\ 
2C(P) 



where C(P) is the number of cells that are predicted to commit by our critical index for 
parameters P; and AP is a small perturbation of all the parameters. 
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